arXiv:1508.00955v4 [astro-ph.SR] 14 Apr 2016 


A Large Convective-Core Overshoot in Kepler 
Target KIC 11081729 

Wuming Yang' 

' Department of Astronomy, Beijing Normai University, Beijing 100875, China 
*E-mail: yangwuming@bnu.edu.cn, yangwuming@ynao.ac.cn 
Received ; Accepted 

Abstract 

The frequency ratios rpi and no of KIC 11081729 decrease firstiy and then increase with the 
increase in frequency. For different spectroscopic constraints, aii modeis with overshooting 
parameter 6ov less than 1.7 can not reproduce the distributions of the ratios. However, the 
distributions of the ratios can be directiy reproduced by modeis with (5ov in the range of about 
1.7 -1.8. The estimations of mass and age of the star can be affected by spectroscopic 
resuits, but the determination of the Sov is not dependent on the spectroscopic resuits. A iarge 
overshooting of convective core may exist in KIC 11081729. The characteristics of roi and 
rio of KIC 11081729 may result from the effects of the large overshooting of convective core. 
The distributions of m and rio of different stars with a convective core can be reproduced 
by the function B{vn,i)- If the value of the critical frequency vo is larger than the value of 
frequency of maximum oscillation power Vmax, a star may have a small convective core and 
(5ov. But if the value of uq is less than that of Umax, the star may have a large convective core 
and i5ov. The function aids in determining the presence of convective core and the size of the 
convective core including overshooting region from observed frequencies. The determination 
is not dependent on the calculation of stellar models. 
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1 Introduction 


Asteroseismology has proved to be a powerful tool for determining the fundamental parameters 
of stars, diagnosing the internal structure of stars, and probing physical processes in stellar interi¬ 
ors by comparing observed oscillation characteristics with those calculated from theoretical models 
(Mazumdar et al. 2006; Cunha & Metcalfe 2007; Christensen-Dalsgaard & Houdek 2010; Yang et al. 
2012; Silva Aguirre et al. 2013; Liu et al. 2014; Guenther et al. 2014; Chaplin et al. 2014; Metcalfe 
et al. 2014; Tian et al. 2014). 

However, the stellar model found solely by matching individual frequencies of oscillation 
does not always properly reproduce other characteristics of oscillations (Deheuvels et al. 2010; Silva 
Aguirre et al. 2013; Liu et al. 2014), such as the ratios of small-to-large separations, rio and roi, which 
are sensitive to conditions in stellar core (Roxburgh & Vorontsov 2003). The small separations are 
defined as (Roxburgh & Vorontsov 2003) 

dlo{n) = + 2l/n,l - l^n+l,o) (1) 

and 

<^01 (^) = ~ ( 2 ) 

In calculation, equations (1) and (2) are rewritten as the smoother five-point separations. 

Stars with a mass larger than 1.1 Mq may have a convective core during their main sequence 
(MS). The overshooting of the convective core extends the region of material mixing by a distance 5ov 
Hp above the top of the convective core that is determined by Schwarzchild criterion, where Hp is the 
local pressure scale-height and 5ov is a free parameter. The overshooting brings more H-rich material 
into the core, prolongs the lifetime of the burning of core hydrogen, and changes the internal structure 
of stars and the global characteristics of the following giant stages (Schroder et al. 1997; Yang et al. 
2012). The determination of the global parameters of stars by asteroseismology or other studies based 
on stellar evolution can be directly affected by the overshooting (Mazumdar et al. 2006; Yang et al. 
2012 ). 

However, for stars with a mass of around 1.1 Mq, there may or may not exist a convective 
core in their interior, depending on the input physics used in the computation of their evolutions 
(Christensen-Dalsgaard & Houdek 2010). Up to now, there is no a direct method to determine the 
presence and the size of the convective core and the extension of overshooting from observed data. 
Therefore, finding a method of determining the presence and the size of the convective core is impor¬ 
tant for understanding the structure and evolution of stars. 

Deheuvels et al. (2010) argued that the presence or absence of a convective core in stars can 
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be indicated by the small separations and that the star HD 203608 with M ~ 0.94 Mq and t ~ 6.7 
Gyr in fact has a convective core. De Meulenaer et al. (2010) studied the oscillations of a Centauri 
A (1.105 ± 0.007 Mq) and concluded that the doi allows them to set an upper limit to the amount of 
convective-core overshooting and that the model of a Centauri A with a radiative core reproduces the 
observed roi significantly better than the model with a convective core. Moreover, Silva Aguirre et al. 
(2013) tried to detect the convective core of KIC 6106415 and KIC 12009504. They obtained that the 
mass of KIC 6106415 is 1.11 ±0.05 Mq and that of KIC 12009504 is 1.15 ±0.04 Mq, and concluded 
that a convective core and core overshooting exist in KIC 12009504, but could not determine whether 
a convective core exists in KIC 6I064I5. Tian et al. (2014) studied the oscillations of KIC 6225718 
and concluded that either a small convective core or no convective core exists in the star. 

For MS stars, the ratios generally decrease with frequency. Liu et al. (2014) argued that the ra¬ 
tios affected by the overshooting of convection core could exhibit an increase behavior. By using this 
seismic tool, the value of overshooting parameter (5ov is restricted between 0.4 and 0.8 for HD 49933 
(Liu et al. 2014) and between 1.2 and 1.6 for KIC 2837475 (Yang et al. 2015). Moreover, Guenther 
et al. (2014) found that the value of 5ov is between 0.9 and 1.5 for Procyon, which is supported by 
the work of Bond et al. (2015). The values of (5ov are larger than the value of 0.1 — 0.3 estimated 
by comparing the theoretical and observational color-magnitude diagram of star clusters (Prather & 
Demarque 1974; Demarque et al. 1994) and characteristics of eclipsing binary stars (Schroder et al. 
1997). 

However, the presence of the large overshooting is supported by the convective theory of 
stars and numerical simulation. For example, Xiong (1985) shows that if the overshooting distance 
is defined as the distance up to which mixing of matter extends, the overshooting distance of the 
convective core can reach 1.4ci Hp in his convective theory, where the value of the ci is between 1/3 
and 1; the numerical simulation of Tian et al. (2009) shows that the penetration distance of downward 
overshooting of convection of giant stars can reach 1 — 2 Hp. Moreover, Valle et al. (2016) show that 
the calibration of the convective core overshooting of double-lined eclipsing binaries with the mass in 
the range of 1.1 — 1.6 Mq is poorly reliable. They suggested that asteroseismic data may be required 
for the calibration of 6ov 

KIC 11081729 is an F5 star (Wright et al. 2003). The frequencies of p-modes of KIC 
11081729 have been extracted by Appourchaux et al. (2012). Combing asteroseismical and non- 
asteroseismical data, several authors (Huber et al. 2014; Metcalfe et al. 2014) have investigated the 
fundamental parameters of KIC 11081729. The mass determined by Huber et al. (2014) is 1.36lo gg 
Mq for KIC 11081729, but that determined by Metcalfe et al. (2014) is 1.26 ±0.03 Mq. The observed 
ratios rio and rgi of KIC 11081729 decrease firstly and then increase with the increase in frequencies. 
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The characteristics may derive from the effects of overshooting of convective core. If these character¬ 
istics can be directly reproduced by stellar models, which may aid in understanding the origin of the 
characteristics and confirming the presence of the large overshooting of convective core. 

In this work, we focus mainly on examining whether the observed characteristics of KIC 
11081729 can be directly reproduced by stellar models and can be explained by overshooting of 
convective core. In order to find the best model of KIC 11081729, the chi squared method was used. 
First, an approximate set of models was found out with the constraints of luminosity, Teff, and [Fe/H]. 
Then around this set of solutions, we sought for best models that match both non-seismic constraints 
and the individual frequencies extracted by Appourchaux et al. (2012), and then we compared the 
observed rio and roi with those calculated from models. And finally, a function i) was deduced 
and tested. In Section 2, stellar models of KIC 11081729 are introduced. In Section 3, the function 
is deduced, and in Section 4, the results are summarized and discussed. 

2 STELLAR MODELS 

2.1 Input Physics 

In order to obtain the best model of KIC 11081729, we construct a grid of evolutionary models 
using the Yale Rotation Evolution Code (YREC) (Pinsonneault et al. 1989; Yang & Bi 2007). For 
the microphysics, the OPAL equation-of-state table EOS2005 (Rogers & Nayfonov 2002), OPAL 
opacity table GN93 (Iglesias & Rogers 1996), and the low-temperature opacity tables of Alexander 
& Ferguson (1994) are used. The standard mixing-length theory is adopted to treat convection. The 
value of the mixing-length parameter a calibrated to the Sun is 1.74 for the YREC. But it is a free 
parameter in this work. The full mixing of material is assumed in overshooting region of models. The 
diffusion and settling of both helium and heavy elements are taken into account in models with a mass 
less than 1.30 Mq by using the diffusion coefficients of Thoul et al. (1994). In the model calculation, 
the initial helium mass fraction is fixed at 0.248 and 0.295. The value of 0.248 is the standard big 
bang nucleosynthesis value (Spergel et al. 2007). The values of input parameters, mass (M), a and 
6ov, and heavy-element abundance {Z) of Zero-Age MS models are listed in Table 1. The models are 
constructed from ZAMS to the end of MS or the subgiant stage. Some of evolutionary tracks of the 
models are shown in Figure 1 as an example. 

The low-/ p-mode frequencies of each model are calculated by using the Guenther’s adiabatic 
oscillation code (Guenther 1994). For the modes with a given degree /, the frequencies i^corrin) 
corrected from the near-surface effects of the model are computed by using equation (Kjeldsen et al. 
2008) 
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Table 1. Input parameters for model tracks. 


The symbol S Indicates the resolution of the parameters. 




Variable Minimum Maximum 

<5 


M/Mq 

1.01 

1.50 

<0.02 

a 

1.65 

2.15 

0.1 

Sov 

0.0 

1.8“ 

0.2 

Zi 

0.010 

0.040 

0.002 


“ Supplemented by ^ov = 1.7. 



Fig. 1. Hertzsprung-Russell diagram of models. Panel (a) shows the evolutionary tracks of models with (5ov = 0 and 0.2 with the spectroscopic constraints of 
Bruntt et al. (2012), while panel (b) depicts those of models with iSov = 1 -8. The circles indicate the most likely models for KIC 11081729. 


f^corr(n) = Umod{n) + (3) 

^max 

where b is fixed to be 4.9, a is a parameter, Uraod{n) is the adiabatie oscillation frequencies of the 
model, i^max is the frequency of maximum power. The value of z/max is 1820 /iHz. The value of a 
is determined from the observed and theoretical frequencies of the modes with the degree I by using 
equations (6) and (10) of Kjeldsen et al. (2008). 

2.2 Spectroscopic Constraints 

The effective temperature of KIC 11081729 may be 6440 K (Wright et al. 2003), 6570^75 K (Ammons 
et al. 2006), 6359 K (Pinsonneault et al. 2012), 6630 ± 70 K (Bruntt et al. 2012), or 6400 ± 127 K 
(Molenda-Zakowicz et al. 2013). The value of [Fe/H] for KIC 11081729 is O.lO^g Js (Ammons et al. 
2006). Combining the value of {Z/X)q = 0.023 of the Sun given by Grevesse & Sauval (1998), the 
ratio of surface metal abundance to hydrogen abundance, {Z/X)^, is estimated to be approximately 
between 0.024 and 0.047. The value of [Fe/H] determined by Bruntt et al. (2012) is —0.12 ± 0.06, 
which corresponds to a {Z/X)^ approximately between 0.015 and 0.020; while that estimated by 
Molenda-Zakowicz et al. (2013) is —0.19 ± 0.22, i.e., the value of {Z/X)s is approximately between 
0.009 and 0.025. Molenda-Zakowicz et al. (2013) have shown that the estimated atmospheric parame¬ 
ters of stars hotter than 6,000 K are affected by the method of spectral analysis. Thus the atmospheric 
parameters determined by Ammons et al. (2006), Bruntt et al. (2012), and Molenda-Zakowicz et al. 
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(2013) were considered in this work. 

The visual magnitude of KIC 11081729 is 9.048 ±0.054 (Droege et al. 2006; Ammons et 
al. 2006). The value of bolometric correction and extinction is estimated from the tables of Flower 
(1996) and Ammons et al. (2006), respectively. The distance of KIC 11081729 is D = pc 
(Ammons et al. 2006; Pickles & Depagne 2010). Then the luminosity of KIC 11081729 is estimated 
to be 4.00 ± 3.28 Lq. The large uncertainty of the luminosity derives from the uncertainty of the 
distance. 

To find the models that can reproduce the observed characteristics of KIC 11081729, we cal¬ 
culated the value of xl of models. The function xl is defined as 


o i=i 


.c: 


'theo 


c: 


obs 




'obs^ 


(4) 


where the symbol C = [L/L0,Teff, {Z/X)^, the presents the observed values of these parame¬ 
ters, while represents the theoretical values. The observational error is given by (t(C°^®). In the 
first step, the models with Xc — f chosen as candidate models. 


2.3 Asteroseismic Constraints 

In order to find the models that can reproduce the oscillation characteristics of KIC 11081729, we 
computed the value of xl for each model. The function xl is defined as 

1 N theo _ ,,obs 

y2 = (5) 

where the quantity u corresponds to the individual frequencies of modes from observation and mod¬ 
elling, is the observational error of the N presents the number of observed modes. The 

value of of corrected frequencies is also computed. The models with xl < 1-0 and xl or xl^^^, 
< 10.0 are chosen as candidate models. 

When a model evolves to the vicinity of the error-box of L/Lq and T^s in the H-R diagram, 
the evolutionary time step of the model is set as small as possible to ensure that consecutive models 
have an approximately equal xl^^^^- some cases, there may be a difference of a few Myr between 
the age of the model with the minimum x^^orr model with the minimum xl- Yang et al. 

(2015) showed that the difference between the model with the minimum xl and the model with the 
minimum xl is insignificant. For a given mass and 5ov, the model with the minimum xl^^rr chosen 
as the best model. 
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Table 2. Parameters of the models with the effective temperture and [Fe/H] determined by Bruntt et 

ai. (2012). The symbol Xc represents the centrai hydrogen abundance of models, while in¬ 
dicates the adiabatic oscillation frequency at which the roi of models reaches a minimum. 


Model 

M 

(Mo) 

TeS 

(K) 

L 

(Lq) 

R 

(Rq) 

age 

(Gyr) 


{Z/X)s 

a 

5ov 

Xc 

xl 

A. l^corr 

Xc 

(mHz) 

Mbl 

1.21 

6524 

3.19 

1.400 

2.475 

0.018 

0.018 

2.05 

0.2 

0.410 

5.5 

5.6 

0.8 

2255 

Mb2 

1.22 

6524 

3.20 

1.403 

2.488 

0.018 

0.017 

2.05 

0.2 

0.417 

5.9 

5.5 

0.8 

2255 

Mb3 

1.23 

6534 

3.22 

1.403 

2.222 

0.018 

0.016 

1.95 

0.2 

0.450 

5.9 

5.6 

0.9 

2346 

Mb4 

1.24 

6556 

3.27 

1.403 

2.073 

0.018 

0.015 

1.95 

0.2 

0.461 

6.6 

5.5 

0.7 

2347 

Mb5 

1.25 

6596 

3.35 

1.403 

1.856 

0.018 

0.014 

1.95 

0.2 

0.479 

7.1 

5.3 

0.6 

2347 

Mb6 

1.26 

6641 

3.47 

1.409 

1.835 

0.018 

0.015 

2.05 

0.2 

0.476 

6.1 

5.4 

0.4 

2347 

Mb? 

1.27 

6555 

3.33 

1.416 

1.907 

0.020 

0.018 

1.95 

0.2 

0.478 

5.4 

5.2 

0.4 

2347 

Mbs 

1.28 

6583 

3.39 

1.416 

1.750 

0.020 

0.017 

1.95 

0.2 

0.492 

7.1 

5.3 

0.2 

2347 

Mb9 

1.29 

6599 

3.46 

1.426 

1.697 

0.020 

0.019 

2.05 

0.0 

0.389 

4.9 

5.4 

0.1 

3562 

MblO 

1.30 

6553 

3.40 

1.432 

2.452 

0.016 

0.022 

1.95 

0.0 

0.327 

5.5 

5.5 

0.9 

4010 


2.3.1 The models with the effective temperature and [Fe/H] of Bruntt 

With the constraints of the effective temperature and the [Fe/H] of Bruntt et al. (2012), we calculated 
the values of Xc> xl^ xl^^„ of o^oh model. Table 2 lists the models that minimize xl + X^^orr f°^ ^ 
given mass. The value of 5ov of the models is less than 0.4. The model Mb8 has the minimum xl + 
xl , which seems to hint that Mb8 is the best model. 

Using the posterior probability distribution function, the mass estimated from the sample of 
the models with xl < 1-0 and X^^orr ^ f^-O is 1.26 ± 0.03 Mq for KIC 11081729, where the error 
bar indicates the 68% level confidence interval. With the same constraints, Metcalfe et al. (2014) also 
obtained the same mass, i.e. 1.26 ±0.03 Mq. The radius is estimated to be 1.41 ±0.01 Rq that is also 
close to 1.382 ± 0.021 Rq of Metcalfe et al. (2014). However, the age of 1.9 ± 0.4 Gyr is larger than 
0.86 ± 0.21 Gyr given by Metcalfe et al. (2014). This is due to the fact that the effects of the settling 
of heavy elements and overshooting of convective core were considered that were not included in the 
models of Metcalfe et al. (2014). This indicates that the estimated age of stars can be significantly 
affected by physical effects that are considered in models. 

However, the calculations show that the models with (5ov <0.2 can not reproduce the distri¬ 
butions of the observed rio and rpi of KIC 11081729 (see Figure 2). This shows that the internal 
structures of these models do not match that of KIC 11081729. The observed ratios change dramat- 
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Fig. 2. The distributions of the observed and theoretical ratios rio and roi as a function of frequency. The theoretical ratios are computed from the corrected 

osciliation frequencies i/oorr of the modeis listed in Table 2. 
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Fig. 3. The distributions of ratios rio and roi of the models listed in Table 3. 


ically at high frequencies, which may derive from the fact that the observed frequencies are affected 
by large-line widths at high frequencies in the same way as that of KIC 6106415 and KIC 12009504 
(Silva Aguirre et al. 2013). Silva Aguirre et al. (2013) suggested that the ratios at high frequencies 
should be excluded from the set of constraints. 

Liu et al. (2014) argued that effects of overshooting of convective core can lead to the fact that 
ratios rio and rgi exhibit an increase behavior. In order to test whether the characteristics of rio and 
roi of KIC 11081729 can be reproduced by the effects of overshooting, we computed the evolutions 
of models with 5ov as large as 1.8. The calculation shows that models with (5ov < 1-7 can not reproduce 
the distributions of the observed rio and roi. Table 3 lists the models that minimize xl + X^<.orr ^ 
given mass and (5ov Model Mb 19 has the minimum xl + X^^orr- Moreover, Figure 3 shows that the 
observed ratios can be reproduced well by models with 5ov = 1-8 and mass in the range of 1.21 — 1.27 
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Table 3. Parameters of the models with the effective temperture and [Fe/H] determined by Bruntt et 

ai. (2012). The symbol Xc represents the centrai hydrogen abundance of models, while in¬ 
dicates the adiabatic oscillation frequency at which the roi of models reaches a minimum. 


Model 

M 

TeS 

L 

R 

age 


{Z/X), 

a 

5ov 

X, 

xl 

zviZcorr 

Xc 

< 


(Mo) 

(K) 

(Le) 

{Rq) 

(Gyr) 









(mHz) 

Mbit 

1.19 

6560 

3.19 

1.384 

4.415 

0.020 

0.015 

1.95 

1.7 

0.566 

20.8 

8.9 

0.8 

818 

Mbl2 

1.21 

6529 

3.20 

1.400 

4.567 

0.024 

0.020 

2.05 

1.7 

0.557 

10.6 

8.1 

0.9 

818 

Mbl3 

1.23 

6605 

3.36 

1.400 

3.793 

0.022 

0.017 

2.05 

1.7 

0.574 

13.7 

8.0 

0.1 

817 

Mb 14 

1.25 

6596 

3.38 

1.409 

3.603 

0.024 

0.019 

2.05 

1.7 

0.576 

12.4 

7.4 

0.2 

817 

Mbl5 

1.27 

6585 

3.40 

1.419 

3.426 

0.026 

0.022 

2.05 

1.7 

0.579 

9.8 

7.0 

0.7 

817 

Mbl6 

1.29 

6582 

3.40 

1.419 

2.856 

0.026 

0.021 

1.95 

1.7 

0.594 

14.4 

7.5 

0.4 

816 

Mbl7 

1.21 

6542 

3.22 

1.400 

4.631 

0.024 

0.020 

2.05 

1.8 

0.559 

10.8 

7.8 

0.7 

818 

Mbl8 

1.23 

6536 

3.22 

1.403 

4.000 

0.024 

0.019 

1.95 

1.8 

0.574 

14.9 

7.7 

0.6 

817 

Mbl9 

1.25 

6607 

3.40 

1.409 

3.654 

0.024 

0.019 

2.05 

1.8 

0.578 

12.2 

7.3 

0.1 

817 

Mb20 

1.27 

6640 

3.49 

1.413 

3.210 

0.024 

0.018 

2.05 

1.8 

0.588 

14.0 

7.5 

0.1 

816 

Mb21 

1.29 

6637 

3.49 

1.416 

2.655 

0.024 

0.017 

1.95 

1.8 

0.602 

22.1 

00 

bo 

0.1 

815 


Model Mb 19 not only has the minimum + X^^orr when 6ov is in the range of 1.7 —1.8, it also 
reproduces the observed ratios. This shows that the structure of model Mb 19 is similar to that of KIC 
11081729, and KIC 11081729 may have a large overshoot of convective core. The mass, radius, and 
age of KIC 11081729 estimated from models with the large (5ov and with Xc < 1-0 < 10.0 

are 1.24 ± 0.03 Mq, 1.40 ± 0.01 Rq, and 3.6 ± 0.7 Gyr, respectively. The mass and radius are close 
to those estimated by Metcalfe et al. (2014). But the age is larger than that determined by Metcalfe et 
al. (2014), which is due to the fact that overshooting of convective core brings more H-rich material 
into the region of hydrogen burning. 

The value of x^<.orr ^ ^ larger than that of Mb8. Moreover, the effective temperature and 

[Fe/H] determined by Molenda-Zakowicz et al. (2013) and Ammons et al. (2006) are different from 
those estimated by Bruntt et al. (2012). The determination of atmospheric parameters can be affected 
by spectroscopic methods (Molenda-Zakowicz et al. 2013). Thus we calculated the models with the 
effective temperature and [Fe/H] determined by Molenda-Zakowicz et al. (2013) and Ammons et al. 
(2006). 
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Table 4. Parameters of the models with the effective temperture and [Fe/H] determined by Molenda- 
Zakowicz et al. (2013). The symbol Xc represents the central hydrogen abundance of models, while 
indicates the adiabatic oscillation frequency at which the roi of models reaches a minimum. 


Model 

M 

TeS 

L 

R 

age 


{Z/X)s 

a 

Sov 

Xe 

xl 

y2 

/X-Ucorr 

Xc 



(Mq) 

(K) 

(Lq) 

(Rq) 

(Gyr) 









(mHz) 

Mml 

1.17 

6291 

2.74 

1.393 

6.738 

0.030 

0.029 

1.95 

1.7 

0.524 

00 

bo 

8.3 

1.0 

819 

Mm2 

1.19 

6326 

2.82 

1.400 

6.070 

0.030 

0.028 

1.95 

1.7 

0.533 

8.5 

8.0 

0.8 

818 

Mm3 

1.21 

6403 

2.99 

1.406 

5.564 

0.030 

0.029 

2.05 

1.7 

0.537 

8.5 

7.7 

0.7 

818 

Mm4 

1.23 

6476 

3.15 

1.413 

4.707 

0.028 

0.026 

2.05 

1.7 

0.553 

8.3 

7.6 

0.5 

818 

Mm5 

1.25 

6509 

3.24 

1.416 

4.180 

0.028 

0.025 

2.05 

1.7 

0.563 

8.4 

13 

0.6 

817 

Mm6 

1.27 

6543 

3.33 

1.421 

3.706 

0.028 

0.025 

2.05 

1.7 

0.572 

8.4 

6.9 

0.8 

817 

Mm7 

1.29 

6535 

3.35 

1.429 

3.527 

0.030 

0.027 

2.05 

1.7 

0.575 

7.6 

6.7 

0.9 

817 

Mm8 

1.17 

6389 

2.89 

1.390 

6.048 

0.026 

0.023 

1.95 

1.8 

0.539 

10.5 

8.4 

0.2 

818 

Mm9 

1.19 

6423 

2.98 

1.396 

5.421 

0.026 

0.023 

1.95 

1.8 

0.549 

10.6 

8.2 

0.2 

818 

MmlO 

1.21 

6457 

3.08 

1.403 

5.296 

0.028 

0.025 

2.05 

1.8 

0.546 

8.6 

7.9 

0.5 

818 

Mmll 

1.23 

6489 

3.17 

1.409 

4.773 

0.028 

0.025 

2.05 

1.8 

0.556 

8.5 

7.8 

0.5 

818 

Mml2 

1.25 

6441 

3.11 

1.417 

4.388 

0.030 

0.027 

1.95 

1.8 

0.564 

8.8 

7.6 

0.6 

817 

Mml3 

1.27 

6515 

3.26 

1.419 

3.607 

0.028 

0.024 

1.95 

1.8 

0.579 

11.0 

7.1 

0.5 

816 

Mm 14 

1.29 

6466 

3.19 

1.425 

3.256 

0.030 

0.026 

1.85 

1.8 

0.588 

12.1 

7.1 

0.5 

816 



V (mHz) 





p (mHz) 




' (mHz) 


Fig. 4. The distributions of the observed and theoretical ratios rio and roi as a function of frequency. The theoretical ratios are computed from the corrected 

frequencies i/corr of the modeis iisted in Table 4. 
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Fig. 5. The distributions of the observed and theoretical ratios rio and roi as a function of frequency. The theoretical ratios are computed from the corrected 

frequencies i/oorr of the modeis iisted in Table 5. 


2.3.2 The models with the effective temperature and [Fe/H] of Molenda-Zakowicz 

Table 4 lists the models that minimize ^ given mass and 5ov The calculations show 

that the models with 5ov < 1-7 can not reproduce the observed ratios. Thus we do not considere these 
models when we estimate the mass and age of KIC 11081729. When the value of 5ov increases to 1.8, 
the observed ratios can be reproduced by models (see Figure 4). This indicates that the distributions 
of ratios rio and roi are mainly dependent on the effects of overshooting of convective core rather 
than on other effects. 

The mass, radius, and age estimated from the models with 6ov in the range of 1.7 — 1.8 and 
with Xc < 1-0 1.230T0.035 Mq, 1.40±0.01 Rq, and 4.6±0.9 Gyr, respectively. 

The age of these models are larger than that of models with the effective temperature and [Fe/H] of 
Bruntt et al. (2012). Thus spectroscopic results can affect the estimation of age of stars. The values of 
of these models are as large as those of the models with the effective temperature and [Fe/H] of 
Bruntt et al. (2012). Thus these models are not better than the models with the effective temperature 
and [Fe/H] of Bruntt et al. (2012). 


2.3.3 The models with the effective temperature and [Fe/H] of Ammons 

The effective temperature determined by Ammons et al. (2006) is between that estimated by Molenda- 
Zakowicz et al. (2013) and that determined by Bruntt et al. (2012). Table 5 lists the models that 
minimize xl + X^^orr ^ given mass and (5ov The distributions of rio and roi of models with (5ov less 
than 1.7 are also not consistent with those of the observed ratios. When the value of 5ov is in the range 
of 1.7 — 1.8, the observed ratios are reproduced well by the models with the large 6ov (see Figure 5). 

Tables 1—5 show that model Ma6 has the minimum x? + xf, in the calculations. Moreover, 
the observed ratios are reproduced by Ma6. Thus Ma6 is the best-fit model in the calculations. When 
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Table 5. Parameters of the models with the effective temperture and [Fe/H] determined by Ammons et 

al. (2006). The symbol Xc represents the central hydrogen abundance of models, while in¬ 
dicates the adiabatic oscillation frequency at which the roi of models reaches a minimum. 


Model 

M 

TeS 

L 

R 

age 


{Z/X), 

a 

5ov 

X, 

xl 

'XUcorr 

xl 



(Mo) 

(K) 

(Le) 

{Rq) 

(Gyr) 









(mHz) 

Mat 

1.21 

6529 

3.20 

1.400 

4.567 

0.024 

0.020 

2.05 

1.7 

0.557 

10.6 

8.1 

0.7 

818 

Ma2 

1.23 

6476 

3.15 

1.413 

4.707 

0.028 

0.026 

2.05 

1.7 

0.553 

8.3 

7.6 

0.9 

818 

Ma3 

1.25 

6509 

3.24 

1.416 

4.180 

0.028 

0.025 

2.05 

1.7 

0.563 

8.4 

13 

0.5 

817 

Ma4 

1.27 

6543 

3.33 

1.421 

3.706 

0.028 

0.025 

2.05 

1.7 

0.572 

8.4 

6.9 

0.4 

817 

Ma5 

1.29 

6536 

3.35 

1.429 

3.527 

0.030 

0.027 

2.05 

1.7 

0.575 

7.6 

6.7 

0.3 

817 

Ma6 

1.30 

6597 

3.52 

1.439 

4.209 

0.034 

0.051 

2.15 

1.7 

0.554 

5.5 

4.8 

0.7 

814 

Ma7 

1.32 

6630 

3.61 

1.442 

3.735 

0.034 

0.051 

2.15 

1.7 

0.563 

5.1 

4.6 

0.9 

814 

Ma8 

1.34 

6592 

3.54 

1.444 

2.939 

0.034 

0.051 

1.95 

1.7 

0.583 

10.9 

6.4 

0.7 

814 

Ma9 

1.21 

6542 

3.22 

1.398 

4.631 

0.024 

0.020 

2.05 

1.8 

0.559 

10.8 

7.8 

0.7 

818 

Mato 

1.23 

6490 

3.15 

1.406 

4.296 

0.026 

0.022 

1.95 

1.8 

0.568 

11.6 

7.7 

0.9 

817 

Mall 

1.25 

6562 

3.32 

1.413 

3.934 

0.026 

0.022 

2.05 

1.8 

0.572 

10.0 

7.2 

0.5 

817 

Mal2 

1.27 

6515 

3.26 

1.419 

3.607 

0.028 

0.024 

1.95 

1.8 

0.579 

11.0 

l.\ 

0.6 

816 

Mal3 

1.29 

6507 

3.28 

1.426 

3.424 

0.030 

0.026 

1.95 

1.8 

0.582 

9.0 

7.0 

0.5 

816 

Mal4 

1.30 

6504 

3.31 

1.435 

5.298 

0.026 

0.036 

1.95 

1.8 

0.599 

7.4 

5.7 

0.3 

815 

Mal5 

1.32 

6487 

3.32 

1.443 

5.054 

0.028 

0.039 

1.95 

1.8 

0.602 

6.5 

5.8 

0.5 

815 

Mal6 

1.34 

6476 

3.33 

1.451 

4.724 

0.030 

0.042 

1.95 

1.8 

0.604 

6.0 

5.8 

0.7 

814 

Mai lb 

1.25 

6522 

3.26 

1.415 

4.240 

0.028 

0.025 

2.05 

1.8 

0.565 

8.5 

lA 

0.5 

817 


the value of 6ov is equal to 1.8, models Mai4 has the minimum xl + Figure 5 shows that the 

ratios of model Mai4 is almost eonsistent with the observed ones. Therefore, Mai4 is ehosen as the 
best-fit model for 5ov = 1-8. 

The mass, radius, and age estimated from the models with 5ov in the range of 1.7 —1.8 and with 
xl < 1.0 and < 10-0 are 1.270 ± 0.035 Mq, 1.420 ± 0.015 Rq, and 3.8 ± 0.7 Gyr, respeetively. 
The median age of 3.8 Gyr is larger than the ages of models with M = 1.27 Mq listed in Table 5. 
This is related to the faet that the table only lists part of the sample. For example, models Mai 1 and 
Mai lb have the same mass and xl^ und almost the same X^^orr’ Figure 6), however, 

they have different age. 

The ealeulations show that the estimated mass and age rely on speetroseopie results but the 
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Fig. 6. The distributions of ratios r-io and roi of models Mai 1 and Mai 1 b as a function of frequency. 


distributions of rio and tqi are mainly dependent on the effects of overshooting of convective core. 
The estimation of 5ov is not affected by the difference in the spectroscopic results. KIC 11081729 
may have a large overshooting of the convective core. 

3 Fitting equation for the observed and theoretical ratios 

The distributions of rio and roi are dependent on 5ov, which may provide an opportunity to diagnose 
the size of convective core from observed frequencies. 

By making use of the asymptotic formula of frequencies, one can get (Liu et al. 2014) 


rio{i^n,i) ~ 2i^„/iv4oAz/, 


( 6 ) 


where 



( 7 ) 


In this equation, c is the adiabatic sound speed at radius r and R is the fiducial radius of the star; r* 
is the inner turning point of the mode with the frequency Equation (6) indicates that ratio rio is 
dependent on the quantity AqAu that is sensitive to the changes in the adiabatic sound speed c. Figure 
7 shows that the changes can directly affect the quantity AqAu. These hint to us that ratios rio and 


roi of stars with a convective core could be affected by the changes in the sound speed. 


The p-mode oscillations of stars are considered to be acoustic standing waves. The variation 
with time t of a standing wave is proportional to cos (cut), where the ui is angular frequency. The 
variation of the modes that penetrate into overshooting region can be affected by convection. For 
I = 1 modes, there should be a critical angular frequency uq. When the angular frequency of modes 
is larger than ojq, the modes penetrate into overshooting region. We assume that the variation with 
time of the acoustic wave that is affected by convection is related to —v4cos(a;of), where A is a free 
parameter. Due to the fact that oscillation frequencies of stars can be observed, thus we want to know 
that the effects of convection core on the characteristics of oscillations vary with frequency rather than 


13 














Fig. 7. Radial distributions of adiabatic sound speed and quantity AqAi/ of models. 


with time. This can be achieved by making use of Fourier transform. Thus one can obtain 

B{oj) =- JQ^''^‘^°[—Acos{ujot)cos{oot)]dt 


rSin(- 


_1 I ^ 

TT L (jJ+Uq ' 

2Au) f 2 N'ku)'' 


( 8 ) 


7r(ajQ—oj^) ^ ojo ’ 

where N is an integer. Taken TV as 1 and uJua = equation (8) can be rewritten as 


= 


2Auri,i . ,27rz/„,i 

sm(——)+5o, 

^n,l) ^0 


27r^(A 


( 9 ) 


where Bq is a constant. The effects of convection core on the characteristics of oscillations should be 
shown by function B{i/nA- 

The value of B{un,i) reaches maxima at around Si^o/S and 7i7o/4 , but reaches minima at 
around z/q and 9z/o/4. Thus the value of i) decreases with frequency between about 3z/o/8 and 
uq, but increases with frequency between about uq and 7z7o/4. The value of B{unA decreases firstly 
and then increases with frequency in the range of about 7i7o/4 — lli^o/4. These characteristics are 
consistent with those of rio and roi calculated from models with a convective core. Thus the ratios 
rio and roi of stars with a convective core could be described by equation (9). 

The observed ratios of KIC 11081729 decrease firstly and then increase with frequency and 
have a minimum between 1804 and 1846 /iHz. Thus one can estimate the value of z/q is between 
about 800 and 820 /iHz. Due to the fact that the value of sin(27rz/„ i/z/q) is 0 at z / o /2 and 3z/o/2, hence 
Bq is not a free parameter. For a given z^g, the value of Bq takes the value of ratio rig or rgi at z/o/2 
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V (/LtHz) 

Fig. 8. The distributions of ratios rio and roi of modeis and quantity B as a function of frequency. The quantity B is computed by using equation (9) with 

A = SOrr, Bq = 0.023, and = 820 /xHz 


or 3izo/2. Thus we obtained the value of Bq is between 0.019 and 0.026 from the ratios of Ma 6 . In 
addition, using equation (9) and the observed z/n,i and rio, the values of parameters A, Bq, and lyg 
are estimated to be 53 ± 21 vr, 0.018 ± 0.003, and 795 ± 21 /iHz,respectively. Figure 8 shows that the 
ratios of KIC 11081729 are reprodueed well by equation (9) with A = 507r, Bq = 0.023, and z/q = 820 
juHz. The distributions of rio and rgi of Ma 6 and Mal4 are reproduced well by equation (9). The 
ratios rio and roi of models Ma 6 and Mal4 reach the minimum at about 810 juHz that are consistent 
with 795 ± 21 /iHz estimated by using equation (9). 

3.1 Test the equation by other stars 

Tian et al. (2014) studied the oscillations of KIC 6225718, and determined that the mass of KIC 
6225718 is l.lOlg 03 Mq, but could not determine whether a convective core exists in KIC 6225718. 
From observed frequencies, one can obtain A = 745 ±527 vr, Bq = 0.021 ±0.008, and izq = 5700 ±423 
/iHz for KIC 6225718. The panel (a) of Figure 9 shows that the observed ratios of KIC 6225718 are 
reproduced well by equation (9) with A = 745 vr, Bq = 0.021, and z/q = 5700 /iHz. Using equation 
(16) of Liu et al. (2014) with fo = 2, the value of zzq for the best model 14 of KIC 6225718 (Tian et al. 
2014) is estimated to be 5890 /iHz. According to equation (9), ratios rio and roi reach the minimum 
at around z/q. The ratios rio and roi calculated from the model 14 reach the minimum at 5642 /iHz. 
The two values are consistent with 5700 ±423 /iHz. The ratios of KIC 6225718 mainly exhibit a 
decreasing behavior. But the observed ratios have a maximum at about 1926 /iHz, which is consistent 
with that B{i/) has a maximum at around 3z/o/8. The model 14 has a small convective core and 5ov 
= 0.2 (Tian et al. 2014). These indicate that a small convective core may exist in KIC 6225718. The 
value of z/max of KIC 6225718 is 2031 /iHz that is much less than the value of uq (5700 ± 423 /iHz). 
The values of A, Bq, and uq computed from the observed frequencies of HD 49933 (Benomar 
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et al. 2009) are 54 ± 14 vr, 0.040 ± 0.006, and 1920 ±46 yuHz. Using equation (16) of Liu et al. 
(2014) with /o = 2, the value of z/q for model M52 of HD 49933 (Liu et al. 2014) is estimated to be 
about 1900 /iHz. The ratios of model M52 reaeh the minimum at about 2057 /rHz. The panel (b) 
of Figure 9 shows that the ratios of model M52 are eonsistent with those eomputed by equation (9) 
with A = ?)l'n, Bq = 0.028, and vq = 1966 /iHz. The model M52 has a large eonveetive eore and 5ov 
= 0.7 (Liu et al. 2014). The value of z /mav of HD 49933 is 1760 /iHz (Appourehaux et al. 2008) that 
is elose to 1920 ± 46 /iHz of uq. The values of extraeted frequeneies of HD 49933 are in the vieinity 
of z/q. The ratios of HD 49933 deerease and then inerease with frequeney, whieh is eonsistent with 
that ealeulated from equation (9). 

Moreover, the panel (e) of Figure 9 shows that the observed and theoretieal ratios of KIC 
2837475 are reprodueed by equation (9) with A = 43 vr, i?o = 0.033, and z/q = 913 /iHz in the range of 
about 1000 — 2000 /iHz. The values of z/q eomputed from the observed frequeneies of KIC 2837475 is 
913 ± 36 /iHz. The ratios rio and roi ealeulated from model Mal4 of KIC 2837475 (Yang et al. 2015) 
reaeh the minimum at about 850 /iHz. The model has a large eonveetive eore and 5 oy = 1.4 (Yang et 
al. 2015). The value of z /mav of KIC 2837475 is 1522 /iHz that is mueh larger than 913 ± 36 /iHz of 
uq. The values of extraeted frequeneies of KIC 2837475 are in the range between about z/q and 9zio/4. 
The ratios of KIC 2837475 inerease and then deerease with frequeney, whieh is also eonsistent with 
that ealeulated from equation (9). 

The distributions of rio and roi of different stars with a eonveetive eore ean be reprodueed 
well by equation (9). The value of z /mav of KIC 6225718 is less than the value of z/q. KIC 6225718 
may have a small eonveetive eore and a small (5ov The value of z/max of HD 49933 is elose to the value 
of z/q. HD 49933 may have a medium (5ov The value of z/max of KIC 2837475 and KIC 11081729 
is larger than that of z/q, respeetively. KIC 2837475 and KIC 11081729 may have a large (5ov The 
equation (9) eould be used to determine the presenee of eonveetive eore and estimate the size of 
overshoot from observed frequeneies. 


16 




















4 DISCUSSION AND SUMMARY 


4.1 Discussion 

Mazumdar et al. (2014) showed that ratios rio and roi can be affected by the glitches at the BCZ and 
the layers of the HellZ. The changes in frequencies caused by the glitch at the BCZ or the layers of 
the HellZ have a periodicity of twice the acoustic depth of the corresponding glitch (Mazumdar et 
al. 2014). The value of the acoustic depth thcIiz and tbcz is about 720 and 2060 s for model Mb8, 
respectively. The value of thsiiz and tbcz of model Mb 19 is 728 and 2100 s, respectively. The 
values of thciiz and tbcz of model Mb8 are approximately equal to those of model Mbl9. Thus 
the glitches should result in a similar effect on the ratios of models Mb8 and Mb 19. However, the 
distributions of ratios rio and roi of model Mb8 are obviously different from those of model Mb 19. 
Thus, the difference between the ratios of Mb8 and those of Mb 19 (see Figures 2 and 3) does not 
derive from the effects of the glitch at the BCZ or the HellZ. 

For different spectroscopic results, only the models with 5ov in the range of about 1.7 — 1.8 
can reproduce the observed ratios. This indicates that the uncertainties of the effective temperature 
and [Fe/H] can not significantly affect the estimation of (5ov This is due to the fact that the estimation 
of 5ov is mainly dependent on the distributions of ratios rio nnd rgi that are only determined by the 
interior structure (Roxburgh & Vorontsov 2003; Oti Floranes et al. 2005). Tables 2 — 5 list the value of 
adiabatic oscillation frequency at which rio nnd roi computed from models arrive at the minimum. 
The values of z/™ of models with 5ov in the range of 1.7 — 1.8 are consistent with that determined by 
equation (9). The calculations show that the value of z/™ is related to the radius of overshooting region 
of models. The larger the radius of the core including the overshooting region, the smaller the value 
of z/™. For KIC 6225718, HD 49933, KIC 2837475, and KIC 11081729 the values of uq estimated 
from observed ratios by using equation (9) are consistent with those obtained from rio nnd roi of the 
best model of the stars. Thus equation (9) aids in understanding the size of the core from the observed 
ratios. If the value of Umax of a star is less than the value of uq estimated by using equation (9), the 
star may have a relatively small convective core. However, if the value of Umax is larger than that of 
uo, the star may have a large convective core. 

The value of 1.7 — 1.8 of the 5ov is much larger than the value adopted usually. However, 
Xiong (1985) shows that the overshooting of convective core could mix material in the distance of 
lAHp. Moreover, the value is consistent with the numerical result for the downward overshooting 
of giants (Tian et al. 2009). A large 5ov means that an efficient mixing takes place in stellar interior. 
Rotation can lead to an increase in convective core, which depends on the efficiency of rotational 
mixing and rotation rate (Maeder 1987; Yang et al. 2013a; Yang et al. 2013b). The effects of rotation 
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explain the extended MS turnoffs of intermediate-age star clusters in the Large Magellanic Cloud, 
where rotational mixing plays an important role (Yang et al. 2013b). The effects of rotational mixing 
and convection overshoot can reconcile the low-Z solar models with helioseismology (Yang 2016). 
Efficient mixing and convection overshoot also may exist in the Sun (Yang 2016). The effect of 
overshooting can mimic the effects of rotational mixing to a certain degree. The surface rotation 
period of KIC 11081729 is about 2.78 days (McQuillan et al. 2014), which is obviously lower than 
the approximately 27 days of the Sun. The large Sov might be related to rotation. 

4.2 Summary 

The observed ratios rio and roi of KIC 11081729 decrease firstly and then increase with frequency. 
There are different spectroscopic results for KIC 11081729, which can affect the determinations of 
mass and age of KIC 11081729 but can not affect the estimation of (iov This is due to the fact that 
the distributions of rio and rgi are mainly dependent on the interior structure. The structure can 
be directly changed by overshooting of convective core. Thus the distributions of rio and roi are 
sensitive to Sqv For the different spectroscopic constraints, the models with Sov less than 1.7 can not 
reproduce the observed ratios. However, the distributions of the observed ratios can be reproduced 
well by models with Sov in the range of about 1.7 — 1.8. A large overshooting of convective core may 
exist in KIC 11081729. The behavior of the ratios of KIC 11081729 may derive from the effects of 
the large overshooting of convective core. 

The ratios rio and roi calculated from models with a convective core reach a minimum at 
about z/q, arrive at a maximum at around 7iyo/4, and then reach a minimum at about 9z/o/4. These 
characteristics can be completely reproduced by equation (9). The distributions of observed and 
theoretical ratios of stars with a convective core can be reproduced well by the equation. The value 
of lyQ decreases with the increase in the radius of convective core and can be estimated by formula 
(9) from observed frequencies. Thus the equation aids in determining the presence and the size of 
convective core including overshooting region from observed frequencies. If the value of of a 
star is less than the value of z/q determined from observed frequencies using equation (9), the star may 
have a small convective core and a small 5ov However, if the value of ly^ax is larger than that of z/q, 
the star may have a large overshooting of convective core. 
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